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ABSTRACT 

We discuss the thermal evolution and Bose-Einstcin condensation of ultra-light dark matter parti- 
cles at finite, realistic cosmological temperatures. We find that if these particles decouple from regular 
matter before Standard model particles annihilate, their temperature will be about 0.9 K. This temper- 
ature is substantially lower than the temperature of CMB neutrinos and thus Big Bang Nucleosynthesis 
remains unaffected. In addition the temperature is consistent with WMAP 7-ycar-|-BAO+H0 obser- 
vations without fine-tuning. We focus on particles of mass of m ~ 10"^'^ eV, which have Compton 
wavelengths of galactic scales. Agglomerations of these particles can form stable halos and naturally 
prohibit small scale structure. They avoid over- abundance of dwarf galaxies and may be favored by 
observations of dark matter distributions. We present numerical as well as approximate analytical 
solutions of the Friedmann-Klein-Gordon equations and study the cosmological evolution of this scalar 
field dark matter from the early universe to the era of matter domination. Today, the particles in the 
ground state mimic presureless matter, while the excited state particles are radiation like. 
Subject headings: (cosmology:) theory, (cosmology:) dark matter, (cosmology:) early universe 



1. INTRODUCTION 

Most matter in the universe is non-luminous. The 
observed flatness of the galactic rotation curves indi- 
cates the presence of dark matter halos around galax- 
ies. Observa tions of the cosmic microwave background 
anisotropics (jSpergel et al. I [20071 ) combined with large- 
scale structu re and type la supernova luminosity data 
(|Reiss et al. 111998. 2004. 20051 iPerimutter et a/.lfTOOl 
constrain cosmological parameters finding that visible 
matter contributes only about 4% of the energy den- 
sity of the universe, as opposed to 22% being dark mat- 
ter and 74% dark energy. More recently, a clear sep- 
aration between the center of baryonic matter and the 
tota l center of mass w as observed in the Bullet clus- 
ter owe et al~\\200($) and la ter in other galaxy cluster 
collisions (|Bradac et al. II2008D . These observations rein- 
force the claim that dark matter is indeed composed of 
weakly interacting particles and is not a modification of 
gravity. 

In the past few decades numerous dark matter can- 
didates have been suggested includin g WIMPs, ax- 
ions, and various spin ze ro bosons (jKamionkowski I 
I2007t iBertone et 011 I2005D . Fundamental spin zero 
particles represented by scalar fields play an impor- 
tant role in particle phy sics models (jPeccei fc Ouinni 
119771: [Torres et aH 120001) . These particles could form 
gravitationally stable structures such as boson stars, 
soliton stars, and galactic halos through some type 
of Jeans inst ability mechanism (ISeidel fc Suen I 119911 : 
lUrena-Lopei1l2002t lAlcubierre et al. 1120031 ). The stabil- 
ity and gravitational wave signatures o f compact scalar 
stars have been stud i ed nu n iericallv (iKhlopov et al. \ 
T985I: ISeidel fc SueTTI \Tmi IBalakrishna et al. I 120061 
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and form gravitationally stable structures are supported 
against collapse by Heisenberg's uncertainty principle 
like boson stars. Structure formation on scales smaller 
than the spreading of an individual boson (the Comp- 
ton wavelength of one par t icle) i s forbidden by quan- 
tum mecha.nics fflu et al. 1 120001 : iSahni fc Wimgl 120001 : 
iLee &: Lim 1 120081 ). Halos formed from ultralight scalars 
with Compton wavelength of galactic scales thus do not 
lead to over-abundance of dwarf galaxies unlike cold dark 
matt e r simulations with heavier bos o ns (iNavarro et aL " 
19961: IMatos fc Urena-Lopez I 120011: lAlcubierre et al. 



200a ISalucci et al. II2003D 



Scalar field halos have been fit to r o tation 
curves of spiral galaxies (iSchunck fc Liddlel 119971: 



Guzman & Urcna-Lopcz ' '2003: Matos & Urcna-LopczJ 
2000; .Arbev et al. , 2001; Bochmcr & Harko... 2007|). 
By using the mass as a fr ee para meter to fit ro- 
tation curves lArbev et al. I (j2001l ) have obtained 
m = (0.4 — 1.6) X 10~^^ eV for n on- interacting ultra- 
light bosons. In a followup paper, I Arbev et al. I (|2002[ ) 
performed a non-thermal analysis of ultralight bosonic 

h alos. 

lUrena-Lopez I (|2009[ ) pointed out that scalars field par- 
ticles can Bose-condense at finite temperatures resurrect- 
ing previ ous work on relativistic Bose-E instein conden- 
sation bv IParker fc Zhang I (j 199 11 11993D . A condensate 
is considered relativistic when the temperature of the 
conden sate is significantly larg er than the mass of one 
boson. IParker fc Zhang 1 (I1993D discuss inflationary ex- 
pansion driven by a relativistic Bose-Einstein condensate 
that then evolves into a radiation dominated universe. 

In this paper we perform a thermal analysis of the post- 
inflationary cosmological behavior of scalar field dark 
matter formed from ultra-light bosons. We use th e quan- 
tum field theory formalism of IParker fc Fulling (|1974f ) 
and extend the analysis that iHu I ( 19821 ) used for a de- 
scription of finite temperature effects in the early uni- 
verse. The bosons are described by a complex scalar 
field to provide a conserved charge. 
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We assume the scalar particles decouple after inflation 
in the early universe, after which the field has a simple 
quadratic potential with no interactions. There are no 
constraints on this field from particle physics or precision 
tests of gravity, as there would be for an interacting light 
scalar field. Ultralight particles (to ^ 10^^'^ eV) form a 
pure ground state Bose-Einstein condensate with a high 
critical temperature behaving like cold dark matter to- 
day. Particles in excited states behave like radiation to- 
day and hence contribute to the amount of hot dark mat- 
ter in the universe, which is constrained by cosmological 
observations. Our model assumes that the interaction of 
the scalar field with normal matter turns off in the early 
universe, at a temperature where most Standard Model 
particles have not yet annihilated, which yields a scalar 
field temperature r$ = 0.9 K with no fine-tuning. A 
reasonable choice of T,s> strongly consistent with cosmo- 
logical observations, is < 1.5 K. Since these estimates 
come from states that behave like radiation, they are in- 
dependent of particle mass. 

In ^J2]we review relativistic Bose-condensation at high 
temperatures and discuss the temperature at which these 
scalars could decouple. fJH] follows the cosmological evo- 
lution of the field. We use units where ft = c = 1. 

2. BOSE-EINSTEIN CONDENSATION AT FINITE 
TEMPERATURE 

We consider a system of ultralight (to ~ 10"^^ eV) 
relativistic bosons represented by complex scalar fields. 
The condition for a relativistic condensate is that the 
temperature of the condens ate T ^ to, which is certainly 
true up to the present day (jUrena- Lopez |[200l . 

In the case of a complex field, there is a conserved 
charge, which is required for traditional Bose-Einstein 
condensation (BEC)^. The charge density is defined as 
the excess of particles n over anti-particles n: 

q = n — n. (1) 

For the excite d states, this charge density is 
(IMukhanov II2005D 



(2) 



where g is the number of degrees of freedom of the system 
and the chemical potential fi < m. The maximum q^-x = 
occurs for ^ = m. For ultra-light bosons, the 
excess of bosons over anti-bosons in the excited states 
is small compared to the number density. Any charge 
added to the system when fj, = m must condense to the 
ground state. Wq is large, the g round stat e is populated 
almost exclusively by particles (jMukhanov 2005). 

The critical temperature below which condensation oc- 
curs is found in terms of the charge density of the dar k 
matter particles (jUrena- Lopez \\2(m : IMukhanov l[2005h 



T<Tr_ = \ — . 



(3) 



When T < Tc the majority of the bosons will condense to 
the ground state. In the ground state the particles will 



behave like non-relativistic matter, while the particles in 
excited states will remain highly relativistic. 

Assuming that BEG occurs and that most particles are 
in the ground state, a first approximation to the total 
dark matter density is 



PDA! ~ in + n)m. 
The density of dark matter today p^ja/ 
pI,, « 23%pc, 
where pc ~ 4.19 x lO^^^eV^. Since n 



and 



Tc 



1.7 X lO^^eV^' 



10^^ K, 



(4) 
(5) 

(6) 
(7) 

which corresponds to a very pure condensate today. Note 
that the required charge density is very high implying the 
necessity of a mechanism that would produce such a large 
asymmetry of scalar particles over anti-particles. 

The matter in the excited states is relativistic and con- 
tributes to the density of hot dark matter phdm, com- 
monly parameterized by the effective number of neutrino 
species Neff. The nominal value due to the three known 
neutrino species plus small cosmological corrections is 
Neff = 3.04. The scalar field temperature r$ can be 
determined from 

PHDM = {Neff - 3.04) I^T^ = , (8) 

with = 1.95 K the neutrino temperature today. 
The value given by the WMAP 7-year results com- 
bined with measurements of the Baryon Acoustic Os- 
cillation and Hubble parameter (iKomatsu et al. Il2010f ) 
is Neff = 4.34t°J^ (68% confidence level). The 68% 
confidence level on Neff corresponds to r$ between 1.51 
K and 2.25 K. 

Big Bang Nucleosynthesis measurements constrain 
relic abundances of deuterium and ^He putting a 2- 
sigma bound of Neff < 3.5 (or Neff < 3.3), which 
depends on the uncertain ty in the ^iJe abundance con- 
sidered (jSteigman I I2007D . This bound corresponds to 
T<j, < 1.5 K (T$ < 1.35 K). Larger values of Neff allow a 
greater contribution of unknown particles like our scalar 
field to Phdm- This could either allow to be larger 
or allow other particles to contribute to phdm- 

We now show that 7$ can be less than the observa- 
tional upper limits and substantially less than T^, with- 
out fine-tuning. Our model assumes that the scalar field 
interactions with normal matter turned off completely 
in the very early universe. Subsequently, the heavier 
Standard Model particles annihilated, dumping their en- 
ergy into the photons, neutrinos, electrons, and positrons 
that remained, increasing their temperature relative to 
that of the scalar field. The photons and leptons have 
about gi — 10.75 degrees of freedom, the complex scalar 
field has 2, and when all the Standard Model particles 
are present g^, ^ 100. Following a standard textbook 
(jMukhanov ll2005|) , entropy conservation gives a relation 
for the temperatures 



^ Recently, ISikivie &: Yangi II2009I ) showed that dark matter ax- 
ions can form a BEC as well. 



-I- 2T| = .g,r, 



(9) 



This gives r$ w 0.9 K, well within our constraints. 



3. COSMOLOGICAL EVOLUTION 

The density evolution of the universe must closely fol- 
low the standard ACDM model at times later than nu- 
cleosynthesis, at a w 10"^*^. During radiation domina- 
tion, our scalar field must be a subdominant contribution 
to the density of the universe, and during matter dom- 
ination it must be a replacement for dark matter. As 
shown below, the excited states are a subdominant con- 
tribution to the density. The macroscopically-occupied 
ground state has p oc at early times and must be con- 
strained to be less dense than the density of radiation for 
at least all times after nucleosynthesis. 

We first derive equations in a general background, as 
well as expressions for density and pressure, then spe- 
cialize to the epochs of radiation and matter domination. 
We also present numerical solutions to the equations and 
show that the initial conditions can be adjusted to satisfy 
the cosmological constraints. 

3.1. Evolution Equations 

The line element in an expanding universe can be writ- 
ten as the Friedmann-Robertson- Walker (FRW) metric 



ds^ = -dt^ + a{tf {dx^ + dy"^ + dz^) . 
The Klein- Gordon equation 

□ $ - = 
is derived from a Lagrangian density of the form 
1 



(10) 



c = 



(11) 



The density and pressure can be defined in the usual 
way 

p= i (at$^at$ -I- + m^^^t) ^ (12) 

where Greek indices vary between 1 and 4 and the index 
j varies between 1 and 3. 

Following iHu I ()1982|) . we perform a series of variable 
transformations to expose the conformal properties of 
the scalar field equation. We introduce a conformal time 
coordinate defined by dt = adr. We also make the sub- 
stitution ^ — a$. The metric is conformally static 



ds' 



i(t)2 [-dr^ + dx^ + dy^ + dz^) 



(13) 



We can now rewrite the Klein-Gordon equation using the 
flat space operator □ = —9^-1-9^-1-9^-1-9^ and remember- 
ing that the d'Alembertian □ = {\/—g)~'^d^{y/—gg^'^du) 



1 



= , 



(14) 



where ' denotes the derivative with respect to t. The last 
two terms in Eq. (jl4p break conformal invariance. The 
invariance could be restored, as in iHu I (|1982( ). by setting 
the mass to zero and adding a term proportional to the 
four-dimensional Ricci scalar. However, in this case the 
field would no longer be minimally coupled. We choose 
to treat the terms that break conformal invariance as a 
perturbation. 



The solution to t he scalar field equation can be decom- 
posed into modes (jParker fc Fullin^ll974[ ) 



(15) 



where H.c. denotes the Hermitian conjugate and i}}^ sat- 
isfies 

a 



dT^ 



2 2 

a m 



^'fe =0 



(16) 



The conserved current in mode k can be written as 
1 



(17) 



The canonical commutation relation of the field $ and its 
conjugate momentum 11 leads to the usual commutation 
relations 

[Aj:,Aj:,]^0, [A,,, A^,] J(fc, P), (18) 

when the conserved current is chosen to be Jq^ = 1 for 
particles and Jok — 1 for anti-particles. The opera- 
tor Aj: corresponds to physical particles and the num- 
ber density of particles is defined to be n =< ^^^^r > 

(jParker fc Fullind[l97l . 

The commutation relations are automatically satisfied 
if we take ipk of the form 



Mr) 



(19) 



Each mode is now characterized by its eigenfunction a;^ 
given by 



d'^ujk 



2 ^d{a^H) 



3/ 2 TT\2 ( 

d{a^H) doJk 



(20) 



da 



4 
^k 



da da 
+ {k'^ + a^rT?)ujl = 0. 

To concomitantly solve the Friedmann equation 

^ttG 



-P 



(21) 



we approximate the Hubble parameter in different epochs 
power laws of the form H = a'/a^ = Hoa~"'. Here p = 
Piad + PA + Pm, the radiation energy density p^ad oc a""*, 
the matter density pm ~ p$ cx a~^ is dominated by dark 
matter and the dark energy term p\ ~ constant in a 
ACDM model. Thus, the exponent is rt = 2 during the 
radiation domination era and n = 3/2 during the matter 
domination era. 

3.2. Radiation domination 

In the radiation domination regime (n = 2) , the scalar 
field Eq. reduces to the fiat space wave equation 
with an effective mass that varies with the scale factor 



dT^ 



(22) 



The mass term is the only perturbation from conformal 
invariance (Since a" = 0.) We extend the analysis of iHu I 
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Fig. 1. — The dark matter density pDM (along with Prad) and 
w = Pdm / Pdm a^rt' displayed as a function of scale factor a. 



()1982D to determine the average density in excited states. 
The density in state k is 



Pk 



(23) 



Modes with k ^ am are effectively massless (cjfe « 
+ a?nn?). The total density in the excited states 

is 



1 



Pk 







cyi-p[LOk/{Ta)\ - 1 
l2^^'"' 



(24) 
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where the terms in Eq. (j23p contribute to higher order. 
When T H, the T"' term dominates the density and 
the excited states behave like radiation. 
For the ground state (fc = 0), Eq. ([20|) can be rewritten 

as 

2dx^ 4 V dx . 



y' + x^y^^O, (25) 
where y and x are dimensionless variables defined by 

(26) 



Ormy, a = 



Ho 



where H = i^orfl"^ with Hor « 1.4 x lO"^^ eV. Note 
that X ~ 1 {a 10^®), which corresponds to H ^ m, is 
the transition to matter-like behavior for these particles. 
When a; ^ 1 (or H ^ m), wc can neglect the x^y^ term. 
Now Eq. has an exact solution of 



Co 



l + C^ix-xoY 



(27) 



where Co and xo are constants. When > 1, y is ap- 
proximately constant. When ^ 1 , the typical behavior 
is that y for small x, peaks at x = xo with a height Cq, 



and then falls off as Cn 



(the higher Co , the narrower 



the peak and hence the transition between the constant 
and x~^ behaviors is more abrupt). Initially, y — con- 
stant and the density po and pressure po for the ground 
state are oc a~^. When x xo and ii x :s> l/|Co| then 



y - -0 
two terms 



^x ^. The pressure and the density then have 



H, 



Po^ 



3/2 
Or 



m' 



4a6mi/2Co 



2H, 



3/2 



(28) 



3/2 
Or 



'Or 



ia^m^/^Co 



2H, 



3/2 
Or 



where the a ^ term dominates at early times. The den- 
sity transitions to a cosmological constant with po = —po 

when a H^^_C^^^^}^-^I^ . If > prad then 
p$ cx (jArbev et al. |[2003 ). However, in this regime 
the above derivation becomes invalid. 

3.3. Matter Domination 

In the matter domination regime {n ~ 3/2), Eq. ([1 
becomes 



dr2 



W2 

2a 



2 2 

a m 



V'fc = 0, 



(29) 



where H = i^oma"^/^ with i?o„ w 7.8 x IQ-^'^eV . 
Using Eq. (|19l) we obtain an equation for ojfe: 



2 



•^k 


duJk 


3a 




4 


da 


+ - 


\da J 




r- 


1 


/ am 


^Om 




2^ ^ 


\Hom 



W2 

-"Om 



(30) 



For the ground state (fc = 0), this equation has an exact 
solution 



Ci sin[4ma3/2/(3Hom) + a] + C2' 



(31) 



where Ci, C2 and the phase a are constants with the 
const raint C9 — (7 ? = 1. This solution is in agreement 
with lArbev et al. I (I2002D . When Ci = the solution 
reduces to wq = aw. Solutions with non-zero Ci oscil- 
late around the uJo = solution. Eq. (|3ip can also be 
written in terms of t as 



Wo 



am 



Ci sm{2mt + a) + C2 



(32) 



The oscillations have a period of n/m a few years 
for TO = 10"^'^ eV). The pressure averages to zero on 
cosmological timescales causing the ground state scalar 
field particles to behave like pressureless matter. 

3.4. Numerical Solution 

Wc also solve Eq. numerically including the effect 
of both radiation and matter in the Hubble parameter, 



H 



Hla-^ 



-3 



(33) 
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The numerical solutions are fully specified by the value 
of the field and its first derivative at a given a, as well 
as an overall scaling of the density. In both cases, the 
overall scaling of the density was chosen to match the 
observed cosmological density of cold dark matter. 

Two representative solutions are shown in Fig. [TJ The 
initial conditions in terms of the x and y variables intro- 
duced above are 2/(10-^) « 0.48 and y'ilO'^) = 0.49 for 
solution 1, and 2/(10"^) sa 5 x 10^ and y'{lO^^) « 5. x lO"^ 
for solution 2. Fig. [IJa) shows the ground state density 
as a function of scale factor, along with the radiation 
density piad- Fig. [T}d displays w = pdm/pdm- The den- 
sity at early times decays like and is determined by 
the mass of the scalar particle and the required density in 
the ground state at late times. Solution 1 requires a large 
density at early times which is not compatible with the 
standard ACDM model. Solution 2 has a phase where 
the ground state density is constant. This allows the 
initial density to be much lower and therefore compati- 
ble with the standard model up to the the approximate 
time of nucleosynthesis (a w 10~^°). The final density 
is the same, but w oscillates at late times around the 
pressureless w = solution. 

4. CONCLUSION 

In this paper we assume dark matter is composed 
of ultralight scalar particles (to — 10~^^ eV) with a 
Compton wavelength of galactic scales. Halos formed 
from such particles naturally do not exhibit small scale 
stru cture, avoiding the ov er-abundance of dwarf galax- 
ies. lUrena-Lopezf ()2009f ) showed that for m < lO"^'^ 
eV Bose-condensation always occurs. The condensate 
has the correct cosmological behavior today. When 
H < m, particles in the ground state behave like presure- 
less matter while particles in excited states act as radia- 



tion. When H > to, the dark matter density is initially 
p$ oc and then switches to a cosmological-constant 
behavior (jArbev et al. |[2002l ) . The density of the excited 
states remains radiation-like p$ ex oc until T ^ H 
(a ^ 10~^^) and is sub-dominant to the density of known 
radiation. 

We find that if these particles decouple from regu- 
lar matter before Standard Model particles annihilate, 
their temperature today is T$ sa 0.9 K. This tempera- 
ture is substantially lower than the temperature of the 
CMB and neutrinos, leaving nucleosynthesis unaffected. 
It is consistent with cosmological constraints on the 
amount of hot dark matter in the universe from WMAP 
7 year-|-BAO-|-H0 observations, which yield a one-sigma 
upper limit of r$ < 2.25 K that is independent of parti- 
cle mass. After decoupling the scalar field has no inter- 
actions and hence cannot be detected by particle physics 
experiments or precision tests of gravity. 

A challenging requirement for our model is the large 
particle- antiparticle asymmetry that may arise from sta- 
tistical fluctuations in the early universe or from asym- 
metric reactions before decoupling. Understanding this 
better is the subject of future work. A more obvious 
challenge is the low mass, which seems unnatural from a 
particle physics perspective. However, even lighter scalar 
fields (to ^ 10~^^eV) have been propo sed to explain the 
accelerated expansion of the universe (jRatra &: Peebles"! 
[T9881: ,Coble et a l. 1997). 
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